Delta and Omega electromagnetic form factors in a three-body covariant 

Bethe-Salpeter approach 
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The electromagnetic form factors of the A and 0. baryons are calculated in the framework of 
Poincare-covariant bound-state equations. The quark-quark interaction is truncated to a single 
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I. INTRODUCTION 

The spatial distribution of hadrons' extensive proper- 
ties, such as mass or electric charge, is of especial rele- 
vance in the understanding of low-energy QCD dynam- 
ics, since they probe the details of the quark-quark and 
gluon-quark interactions. 

The electromagnetic properties of the proton have been 
widely studied experimentally, providing a good picture 
of its internal structure. This is not the case, however, 
for the lightest baryonic resonance, the A(1232). Its 
short lifetime makes the study of its properties difficult 
and only the magnetic moments of A ++ [1—10] and A + 
[10, 11] are known, albeit with large errors. An indirect 
estimation of the A + electric quadrupole moment from 
the 7-/V — > A transition quadrupole moment was given in 
[12]. The 17(1672) decays weakly, instead, and this has 
allowed for a precise measurement of its magnetic dipole 
moment [10]. 

For finite values of the photon momentum the only 
information available comes from lattice QCD calcula- 
tions [13-17]. Although constantly improving, these cal- 
culations suffer from the usual problem of not yet being 
able to work at the physical pion mass. Moreover, the 
limit of vanishing photon momentum is unreachable for 
technical reasons. The calculation of the electromagnetic 
properties of the Delta and Omega baryons has also been 
tackled from a number of constituent quark models [18- 
21], chiral quark-soliton model [22], chiral perturbation 
theory [23, 24] and QCD sum rules [25]. 

In this paper, we investigate the electromagnetic prop- 
erties of the Delta and Omega baryon in the framework 
of covariant Bethe-Salpeter equations. In section II we 
introduce the general formalism of Bethe-Salpeter equa- 
tions (BSE) and Dyson-Schwinger equations (DSE). This 
is followed by a presentation of the truncation used in sec- 
tion III. In section IV the results of our calculation are 
discussed. Finally, we conclude in section V. 



II. BARYON BETHE-SALPETER EQUATION 
AND COUPLING TO AN EXTERNAL FIELD 

The evolution of a three-quark system in quantum 
field theory is described through the six-quark Green's 
function (pi,P2,P3) (in momentum space) or, equiv- 
alently, its amputated version the scattering matrix 
T^{pi,p2,pz). This function can be obtained by solving 
a Dyson equation 1 



T = -iK - iKGqT 



or, equivalently, 



T- 1 = iK' 1 - Go 



t = iTK' 



TG 



(1) 

(2) 
(3) 



where Go is the disconnected product of three full quark 
propagators and —iK is the three-quark interaction ker- 
nel. The latter includes three- and two-particle irre- 
ducible interactions 



K 



(4) 



a=l 



with a denoting the spectator quark (see e.g. Fig. 2). 
The full quark propagator S is obtained by solving the 
quark propagator DSE (see Fig. 1) 

S~\p) = S^ip) + Z lf f rD^(p - q)T v gqq (p, q)S(q) , 

Jq 

(5) 

where So is the (renormalized) bare propagator 

S^{p)^Z 2 {ii> + m) , (6) 

Fg qq is the full quark-gluon vertex and D^ v is the full 
gluon propagator and Z\ / and Z2 are renormalization 
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1 For simplicity, we employ a compact matrix notation in which 
discrete/continuous variables arc implicitly summed/integrated 
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constants. In the Landau gauge the gluon propagator 
takes the form 



k 2 ' 



(7) 



with Z(k 2 ) a dressing function to be determined and T^ v 
the transverse projector 



k k v 



(8) 



The bare quark mass m in (6) must be provided as a 
parameter. 



k=p-q 



-1 



+ -k » 




FIG. 1. Diagrammatic representation of the quark Dyson- 
Schwinger equation (5). Blobs represent fully dressed propa- 
gators or vertices. 

When the three-quark system forms a bound state, the 
Green's function develops a pole at P 2 = — M 2 , with 
P the total quark momentum 



P = Pi + P2 + P3 , 



(9) 



and pi the quark momenta, with M the bound-state 
mass. At the bound state pole one defines 



P 2 + M 2 



(10) 



where AT is a normalization factor which, in the case of 
spin- 3 /2 particles is 2M. The function \t is the bound- 
state Bethe-Salpeter amplitude and , F its charge conju- 
gate. They are expressed as tensor products of flavor, 
color and spin parts which describe a baryon in terms 
of its constituent quarks. For spin- 3 /2 baryons the spin 
part is itself a mixed tensor with four Dirac indices and 
one Lorentz index [26-28]. Substituting (10) in (2) or in 
(3), and keeping only the singular terms, we arrive at the 
Bcthc-Salpeter equation for the three-quark bound state 



= -iKG ^ , 



(11) 



(12) 



A systematic procedure to couple an external gauge 
field to the constituents of a three-particle system de- 
scribed by integral equations is the so-called gauging of 
the equations, introduced in [29-32]. It ensures that the 
resulting equations arc gauge invariant and that there is 
no over-counting of diagrams. For our purposes it suffices 



to say that the gauging of equations acts as a derivative 
on the integral equation. That is, (1) becomes 



= - iK<* - iK^G T 
- iKG^T - iKG Q T» 



(13) 



where the supcrindcx fi denotes a gauged function (that 
is, coupled to the gauge field). This equation can be 
rewritten, using (1), as 

T" = (l + iKGo^x 

{-iK^ - iK^GqT - iKG^T) 
= T (iK~ 1 K t *K~ 1 +G%)T . (14) 

To have an expression for one needs to specify the 
interaction kernel. In the next section we shall obtain the 
gauged kernel in the case of the Rainbow-Ladder trunca- 
tion. The gauged quark propagator allows the introduc- 
tion of the proper vertex through the definition 



= ST^S 



(15) 



which, in the case that concerns us in this paper, repre- 
sents the fully dressed quark-photon vertex. 

The bound-state electromagnetic current J M can be 
introduced at the bound-state pole by 



T (3 )^ — AT- Aft — T**- 



(16) 



Substituting in (14) and using (11) and (12) we arrive at 

J" = * f (-iG K»G + G%)* i . (17) 

The electromagnetic form factors are calculated via the 
identification of (17) with the expression of the current 
imposed by symmetry principles (see Appendix A). 

III. RAINBOW-LADDER TRUNCATION 

To solve (11) one needs to specify the interaction ker- 
nel K. An exact expression for this kernel is in gen- 
eral not available and one has to resort to some trunca- 
tion scheme. The simplest consistent scheme is known as 
Rainbow- Ladder (RL) truncation. This truncation pre- 
serves the axial-vector Ward-Takahashi identity, which 
relates the quark-antiquark interaction kernel and the 
quark-gluon vertex in the quark DSE [33, 34] . In the me- 
son sector this identity ensures that pions are the Gold- 
stone bosons of spontaneous chiral symmetry breaking 
[35] . The RL truncation reduces the quark-antiquark ker- 
nel to a single dresscd-gluon exchange. The full quark- 
gluon vertex is projected onto the tree-level Lorentz 
structure 7^ and the non-perturbative dressing is re- 
stricted to depend on the gluon momentum only and has 
to be modelled. It is customary to include this dress- 
ing and the gluon propagator dressing Z(k 2 ) in a single 
effective interaction a(k 2 ). 
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FIG. 2. Diagrammatic representation of the covariant Faddeev equation in the Rainbow-Ladder truncation (18). 



A. Three-quark bound state equations 

Interactions in the baryon sector are not restricted, in 
principle, by the axial-vector Ward-Takahashi identity. 
However, we adopt here for the quark-quark interaction 
kernel the same truncation scheme and neglect the three- 



particle irreducible interactions. The three-body BSE 
(11) in the RL truncation, which is also known as covari- 
ant Faddeev equation (and, correspondingly, the Bethe- 
Salpeter amplitudes ^ are called Faddeev amplitudes), 
reads (see Fig. 2) 



*a/3 7 i(p, q, P) 



(k) Sp>p>>(k2)Si>y>(ka) *q/3" 7 "x(p (1) ,9 (1) ,-P) 

(—k) S 7 ' 7 "(k 3 )Sa'a"(ki) ^ a "/3-/"I 



(18) 



where we have absorbed the —i factor into the definition 
of K, so that it is defined as 

iW'W = -^C Zl ?*jp. T M „(fc) 7 ^ a , 7 ^, , (19) 

with Z2 the quark renormalization constant. We have 
used the generic index X to refer to the bound state (as 
opposed to the first three Greek indices in the Faddeev 
amplitude which denote the valence quarks); for the case 
of a spin- 3 /2 baryon T consists of a Dirac and a Lorentz 
index. In (18), the flavor part of the Faddeev amplitudes 
has been factored out because the interaction kernel is 
flavor-independent and the factor C = —2/3 stems from 
the traces of the color structures. The Faddeev ampli- 
tudes depend on the quark momenta pi, P2 and p 3 , but 
this dependence can be reexpressed in terms of the total 
momentum P and two relative momenta p and q 

p= (1 -C)P3 - C(pi +P2), Pi = -<?- 1 + 



'i 



pi - pi 



P =Pl+P2+P3, 



P2 — Q P ■ 

P2 q 2 2 

p 3 = p + CP , 



(20) 



with £ a free momentum partitioning parameter, which 
we choose £ = 1/3 for numerical convenience. The in- 
ternal quark propagators depend on the internal quark 



momenta fcj = p, — k and fej = pi + k, with k the gluon 
momentum. The internal relative momenta, for each of 
the three terms in the Faddeev equation, are 



P i3) =P, 



= p + k, p {2) =p-k, 
q {1) =q-k/2, q^=q-k/2, q&>=q+k 



(21) 

The quark DSE in the RL truncation reduces to 

(p) = So,*? (p) + f K aa , pp , {k)S a .p, (q) , (22) 

J q 

where now in K, see (19), we have C = 4/3 and k = p—q. 

B. Bound state electromagnetic current and 
quark-photon vertex 

The expression for the current (17) simplifies consid- 
erably in the RL truncation. Since is absent and 
is reduced to the exchange of a neutral particle, the 
photon can only couple to the quark propagator through 
the term S" 1 in (4). Defining V 1 = (SS^Y = we 
obtain 

(S-y 1 = -S^S^S- 1 = -T" , (23) 
where we used (15). Thus (17) becomes (cf. Fig. 3) 



_ — _ v __ 

VTr a iTr _ _ > - VTr — 



■©— o — e- 



j M = x ^ ^ — 9 — * = - S - * * 



perm. o , perm. 



FIG. 3. Diagrammatic representation of the current (17) in the Rainbow-Ladder truncation, see Eq. (24). 
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(^aP-fliP 



{3} {3} 

i ! y i i 



(3} ^{3} J 3 } 



i>0 > 



(24) 



where we defined 



*a/ 7 I = / 
' k 



(fc) Sp'p" {p 2 - k)Sy~ ( » {p 3 + k) ^ a /3"Y'l{p + k,q-k/2, P) 



(25) 



as a result of the first term in the Faddeev equation (18) 
and in a similar fashion we define *^ 2 > and *^ 3 >. The 
injected momentum Q is introduced via the final and 
initial momenta of the interacting quark k 

pf=p K ±^. (26) 

The relative momenta in the respective terms of (24) are, 
using the definitions in (20), 

p#=*±(i-C)§, = 

and since the initial and final states are on-shell, the total 
momenta are constrained to be PI — Pj — — M 2 , with 
M the mass of the bound-state. As is the case for the 
Faddeev equation, the three terms in (24) are formally 
the same when the momentum partitioning parameter is 
chosen to be £ = 1/3. 

The quark-photon vertex T M can naturally be incorpo- 
rated in the framework of covariant bound-state equa- 
tions by calculating it from an inhomogenous Bethc- 



Salpeter equation 

r"(p,Q) = iZ*f 

+ f K q3 (S(k + Q/2)T"(k,Q)S(k-Q/2)) , 

Jk 

(28) 

and using for K qq the RL kernel (19) with C = 4/3 
and for the quark propagator S the solutions of the RL- 
truncated quark DSE (22). We calculate this in the ap- 
propriate moving frame following Rcf. [36]. 



C. Effective interactions 

The appearance of the effective interaction in (19) will 
a priori introduce a model dependence on our results. In 
fact, this is the only model input of the approach. To 
assess how strong is this dependence and to identify the 
possible model-independent features, we use two different 
models for the effective interaction in our calculations. 

The first model we use is known as the Maris- Tandy 
model [37, 38] and has dominated hadron studies within 
Rainbow- Ladder. This dominance is well-earned since 
this ansatz performs very well when it comes to the 
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purely phcnomcnological calculation of ground-state me- 
son and baryon properties. However, this model has no 
clear connection to QCD in the intermediate- and low- 
momentum regime and is, therefore, not entirely satisfac- 
tory to gain understanding of the formation of hadronic 
bound-states in QCD. On the other hand, with the rapid 
improvement in our knowledge of QCD Green's functions 
from both lattice and functional approaches, it is possible 
to define different effective interactions which, presum- 
ably, capture more faithfully some of QCD's features. 
Based on this, an effective interaction has been proposed 
in Ref. [39]. 

Note that the fact that an effective interaction captures 
more features of QCD does not necessarily mean that it 
will perform better phenomenologically. This is because 
the interaction is used within a given truncation scheme 
and, therefore, if one wants to reproduce hadron proper- 
ties the model has to be tuned to account for the effect 
of the missing contributions. In particular, it has been 
shown in [40] that dynamical quark-mass generation is 
accompanied by the appearance of scalar components in 
the quark-gluon vertex. An application of this beyond 
Rainbow-Ladder has been pursued in Rcfs. [41, 42] with 
a non-diagrammatic means provided in Ref. [43]. In ad- 
dition, unquenching effects in the form of a pion back 
coupling to the quark propagator and two-body kernel 
have been investigated in Refs. [44, 45]. However, none 
of these methods have yet been extended to the covari- 
ant three-body problem presented here and so we restrict 
ourselves to Rainbow-Ladder. Since we lose many com- 
ponents of the quark-gluon vertex we therefore construct 
an effective interaction that attempts to mimic their con- 
tribution. 

In this respect, both models described below are de- 
signed to correctly reproduce dynamical chiral-symmetry 
breaking as well as pion properties at the physical u/d 
mass. This means that they account for missing effects 
in the bound-state pscudoscalar meson sector and at this 
quark mass. As a consequence, both interactions have 
similar strength at the intermediate momentum region 
~ 0.5 — 1 GeV (see Fig. 4). These two interactions have 
previously been compared in Ref. [46]. 



1. Maria-Tandy model 

In the Maris- Tandy (MT) model [37, 38] the effective 
running coupling is given by 



a c s{q 2 ) =7T77 7 



27T7 TO (l - e 



ln[e 



l + (l + g 2 /A& OD ) 2 ] 



(29) 



which reproduces the one-loop QCD behavior in the UV 
and features a Gaussian distribution in the intermediate 
momentum region (see Fig. 4) that provides dynamical 
chiral symmetry breaking. The scale A t = 1 GeV is 



introduced for technical reasons and has no impact on 
the results. Therefore, the interaction strength is char- 
acterized by an energy scale A, fixed to A = 0.74 GeV 
to reproduce correctly the pion decay constant from the 
RL-truncated meson-BSE. The dimensionless parameter 
r\ controls the width of the interaction. For the anoma- 
lous dimension we use 7 m = 12/(1 INc — 2 Nf) = 12/25, 
corresponding to Nf = 4 flavors and N c = 3 colors. For 
the QCD scale Aqcd = 0.234 GeV. Many ground-state 
hadron obscrvables have been found to be almost insen- 
sitive to the value of rj around r) = 1.8 (see, e.g. [47-49]). 
This has been used as an argument in favor of the model 
independence of Rainbow-Ladder results. Instead of pur- 
suing this line of research, we prefer to introduce a new, 
non-related model to evaluate the validity of those asser- 
tions. 

Note that in the numerical resolution of the quark DSE 
we employ the Pauli-Villars rcgularization method of the 
integrals, with a mass scale of 200 GeV. Moreover, for this 
model, we fit the quark masses, at the renormalization 
scale fx = 19 GeV, to be 3.7, 85.2, 869 and 3750 MeV for 
the u/d, s, c, and b quarks, respectively. 



2. Soft-divergence model 



The model of Ref. [39], called soft-divergence or SD 
model from here on, is motivated by the desire to account 
for the (7 J 4(l)-anomaly by the Kogut-Susskind mecha- 
nism [50, 51]. The effective coupling is constructed as 
the product of the gluon dressing [52, 53] and a model 
for the non-perturbative behavior of the quark-gluon ver- 









/ V 


MT 

SD 




/ \ x> 
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0.1 1 

q 2 (GeV 2 ) 



100 1,000 



FIG. 4. Comparison of the Maris- Tandy (MT) (29), with n = 
1.8, and the soft-divergence (SD) (30) effective interactions. 
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tex [40], 



a e ff(<7 2 ) =C 



Ik / x -k-1/2 

x \ I y 



1 + xJ \l + y 

«o + 0[/y a; \ ' ( a uv x 
A 



1 + x 



1 + x 



(30) 



The four terms in parentheses are: the IR scaling of the 
gluon propagator; IR scaling of the quark-gluon vertex; 
logarithmic running of the gluon propagator; and the log- 
arithmic running of the quark-gluon vertex. Addition- 
ally, the last two are constructed to interpolate between 
the IR and UV behavior. The remaining terms are de- 
fined as 



A = 



As 



1 + 2, l + (y-l) 2 
where 



o-uv = 7!"7n 



1 



<?7A 



g 2 /A 2 



YM J 



X 

V — 1 / 1S -IR l 

— 2 I A 2 

z — 1 f IV MOM i 



1 

In z z — 1 

(31) 



(32) 
(33) 
(34) 



and a = 8.915/7V C . Here, A YM = 0.71 GeV is the dy- 
namically generated Yang-Mills scale, while Amom — 0.5 
GeV corresponds to the one-loop perturbative running. 
The IR scaling exponent is n = 0.595353, and the one- 
loop anomalous dimensions are related via 1 + 7o = 
-2S = |JVc7m, with lrn = 12/(11JV C - 2N f ). We 
choose Nf = 5 active quark flavors at the renormaliza- 
tion point \i = 19 GeV. The constant C = 0.968 is cho- 
sen such that a e g runs appropriately in the UV. Finally, 
Air = 0.42 GeV, \ s = 6.25, and X B = 21.83 deter- 
mine the IR properties of the quark-gluon vertex and are 
fitted such that the properties of tt, K and p mesons 
are all reasonably well reproduced. The quark masses at 
fi = 19 GeV arc 2.76, 55.3, 688 and 3410 McV for the 
u/d, s, c, and b quarks, respectively 



3. A remark on missing mesonic effects 

The MT and SD model both rely upon the phe- 
nomenology of dynamical chiral symmetry breaking in 
the light quark sector to determine their parameters. 
Therefore effects we might consider to be beyond RL are 
absorbed into the model parameterisation. In particu- 
lar, since these are determined in the light-quark sector 
we implicitly include those contributions due to interac- 
tions at the hadronic level. Here the pion as the lightest 
hadron plays a special role in the dressing of baryons. 
Amongst these contributions, non-perturbative pionic ef- 
fects - also sometimes called pion-cloud effects, see e.g. 
Ref. [54] and references therein - are expected to have 
a sizeable influence on hadron properties like the masses 



or the decay constants. Consequently when fixing the 
model parameters in the light-quark sector, large parts 
of these so-called pion-cloud contributions are "paramc- 
terised" in c/., the discussion in Ref. [55]. According to 
Zweig's rule the meson-cloud around the triple-strange 
D will be mostly constituted of kaons. Due to their 
higher mass as compared to pions, perturbative as well as 
non-perturbative mesonic effects are significantly smaller 
for the ground state properties of the f2 than for the 
ground state properties of the A. However, as we do not 
change the parameters of the model for the Q we expect 
to find larger deviations from experimental values. This 
is because we actually then overestimate the beyond RL 
effects; they look larger despite being actually smaller. 
This should be kept in mind when comparing our results 
to lattice data and experimental observations. 



IV. RESULTS 

We computed the electromagnetic current of the Delta 
and the Omega baryons and extracted the corresponding 
form factors using (A12). As explained in previous sec- 
tions, the interaction parameters and bare quark masses 
arc 1 fitted to reproduce meson properties. In the baryon 
sector, therefore, there are no further parameters to be 
fixed. 

Of the four A(1232) isospin partners, we restrict the 
discussion to the A + since, due to the assumption of 
isospin symmetry, the form factors of the remaining iso- 
partners can be obtained by multiplying with the corre- 
sponding baryon charge. This, in particular, implies that 
all A form factors are identically zero in our approach. 

The solution of the Faddeev equation, (18), and the 
subsequent calculation of the electromagnetic current via 
(24) is a numerically complicated task, chiefly as a con- 
sequence of the expansion of the Faddeev amplitudes in 
128 Lorentz covariants and in a number of Chebyshev 
polynomials for the angular dependence, which entails 
that one must solve for an equal number of coefficients. 
Due to CPU-time and memory limitations, the number 
of quadrature points used in the numerical integrations 
must be kept small. Moreover, the presence of inverse 
powers of Q in the equations for the extraction of the 
form factors (A12) implies that, to obtain reliable results 
at low Q and even finite results in the limit Q — > 0, very 
delicate cancelations among the many terms that con- 
tribute to the current must take place. For these reasons, 
that limit is difficult to reach with our current resources, 
especially for the electric quadrupole, see also Ref. [56], 
and magnetic octupole form factors. 



A. Electric monopole form factor and charge radius 

The calculated electric monopole form factor Geo(Q 2 ) 
for the A + is shown in the upper-left panel of Fig. 5 and 
compared to lattice calculations using dynamical Wilson 




FIG. 5. Electromagnetic form factors for the A + using the Maris-tandy (MT) and the soft-divergence (SD) models. We 
compare with unquenched lattice data (DWF) at three different pion masses [13, 14] 



fcrmions at three different pion masses [13, 14]. The nat- 
ural scale associated to the problem is the Delta mass; 
since MT and SD models, as well as lattice calculations, 
give different values for this mass, we plot the evolution 
of the form factors in terms of the dimcnsionless quantity 
Q 2 /M 2 to remove the scale ambiguity that appears in the 
comparison of results using different approaches/models. 
We stress again that, since we assume isospin symmetry, 
the form factors for the A ++ , A and A~ are obtained 
by multiplying the former by the corresponding charge. 

We see from Fig. 5 that both the MT and the SD mod- 
els show good agreement with lattice calculations. The 
Q 2 -evolution of Geo differs slightly for the two models 
we considered. However, one must bear in mind that 
we are working here with the simplest chiral-symmetry- 
preserving interaction kernel (namely, the RL kernel). 
Since the effective couplings are tailored to reproduce 
meson observables, we consider it sufficient if they re- 
produce baryon properties at the level of a few percent. 
From this point of view, we can say that the behavior of 



Geo(Q 2 ) is qualitatively model independent in our ap- 
proach. 

The charge radius is calculated using the equation 



< r 



EO 



>= - 



6 dG E0 (Q 2 ) 
G E o(0) dQ 2 



(35) 



and the results are shown in Table I for the MT and SD 
models as well as for lattice calculations. As before, we 
can suppress the scale dependence of the charge radius 
by calculating the dimcnsionless quantity < r E0 > M\. 
This quantity shows a better agreement with the lattice 
data than the dimensionful charge radius does, although 
the value for the SD model is significantly larger. 

It is worth mentioning that chiral perturbation theory 
shows that, when the A — > Nir decay channel opens, 
the charge radius changes abruptly to a lower value [24] . 
Since in our calculation we do not provide a mechanism 
for the Delta to decay, it is therefore reasonable that in 
a full calculation this would lead to a lower result for 
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FIG. 6. Electromagnetic form factors for the SI using the Maris-tandy (MT) and the soft-divergence (SD) models. We compare 
with unquenched (DWF) and mixed (Hyb.) lattice data at three different pion masses [15] 



< r 



£0 



>. This effect, nevertheless, would be compen- 



sated partly by the inclusion of mesonic effects. 



Since we asumc isospin symmetry, in our framework 
the fl baryon is identical to the A but evaluated at a 
different current-quark mass. We show the evolution of 
the electromagnetic form factors for the fi~ in Fig. 6. 
The calculation shows good agreement with lattice data 
for both models and, as before, a qualitative agreement 
between them. The electric charge radius is shown in Ta- 
ble II. In this case the calculated charge radius is smaller 
than the lattice values. However, the dimcnsionlcss quan- 
tity < r% Q > Mq shows good agreement between our 
results and lattice. Also, our result for this quantity 
shows little quark-mass dependence, as can be seen by 
comparing the values for the Q and the A; presumably, 
the inclusion of pion-cloud effects, or indeed other fla- 
vor dependent contributions beyond that of Rainbow- 
Ladder, would account for the quark-mass dependence 
of the charge radius. 



B. Magnetic dipole form factor 

As already mentioned above, the magnetic moments 
of the A + and A ++ arc two of the few electromagnetic 
properties of the Delta for which we have experimental 
input. The value at Q 2 = of the magnetic dipole form 
factor Gmi(0) for the A + , which is related to the mag- 
netic moment via the relation 



fJ-A 



2M A 



Gmi(O) , 



(36) 



is given in Table I. We find good agreement between our 
results and the lattice data at different pion masses. The 
value of Gmi(O) for the fi~ is shown in Table II. Here the 
comparison with lattice is less favorable, and we clearly 
underestimate the experimental value which, in this case, 
is very accurately measured. This is a signature of miss- 
ing meson-cloud effects whose relevance is, as discussed 
in the last section, somewhat obscured at the u/d quark 
mass region since the effective interactions are fitted for 
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F-MT F-SD 


DW1 DW2 DW3 


Exp. 


M A (GeV) 


1.22 1.22 


1.395 (18) 1.559 (19) 1.687 (15) 


1.232 (2) 


< r% > (fm 2 ) 


0.50 0.61 


0.373 (21) 0.353 (12) 0.279 (6) 




< r| > Ml 


0.75 0.91 


0.726 (36) 0.858 (25) 0.794 (14) 




Gmi(O) 


2.38 2.77 


2.35 (16) 2.68 (13) 2.589 (78) 


3.54±2;» 



TABLE I. Comparison of results for the A + mass, charge radius < r% > and for Gmi(0) (oc /j,). We compare our results for the 
MT model (F-MT) and for the SD model (F-SD) with a lattice calculation with dynamical Wilson fermions at = 384 MeV 
(DW1), m n = 509 MeV (DW2) and = 691 MeV (DW3) [13, 14]. For G M i(0) we also compare with the experimental value 
[10, 11]. 





F-MT F-SD 


DW1 DW2 Hyb. 


Exp. 


Mn(GeV) 


1.65 1.80 


1.76 (2) 1.77 (3) 1.78 (3) 


1.672 


< r% > (fm 2 ) 


0.27 0.27 


0.355 (14) 0.353 (8) 0.338 (9) 




< r% > Mq 


0.74 0.89 


0.726 (36) 0.858 (25) 0.794 (14) 




Gmi(0) 


-2.41 -2.71 


-3.443 (173) -3.601 (109) -3.368 (80) 


-3.52 (9) 



TABLE II. Comparison of results for the Q mass, charge radius < r% > and for Gmi(0) (oc fi). We compare our results for the 
MT model (F-MT) and for the SD model (F-SD) with a lattice calculation with dynamical Wilson fermions at m-n- = 297 MeV 
(DW1), m w = 330 MeV (DW2) and with a hybrid action at m n = 353 MeV (Hyb) [15]. For Gmi(0) we also compare to the 
experimental value [10]. 



that sector, thus in a sense incorporating pion-cloud ef- 
fects in the parameters of the model. 

The evolution of Gmi with the photon momentum also 
compares favorably with lattice results in the case of the 
A. Again, this is not the case for the f2 as now both mod- 
els differ significantly from lattice calculations at low-Q 2 , 
where pion- and kaon-cloud effects arc expected to be 
more relevant. 



C. Electric quadrupole form factor 

A non-vanishing value for the electric quadrupole 
moment signals the deformation of the electric charge 
distribution from sphericity. It would be identically 
zero if the baryon were formed only by s-wave compo- 
nents. In our approach, the presence of higher angular- 
momentum components is a natural consequence of re- 
quiring Poincare covariance [27]. Nevertheless, the rela- 
tive importance of these components is dictated by the 
dynamics and we could still obtain a non-trivial vanishing 
value for this moment. 

We show our calculations for the electric quadrupole 
form factor and its evolution with Q 2 in the bottom- 
right panel of Fig. 5. Although the precise vale of 
Ge2(0) is very sensitive to numerical accuracy (due to 
the presence of an 1/Q 4 factor when extracting the form 
factor from the electromagnetic current; see (A12)), we 
clearly see that for both the MT and SD models it is 
non- vanishing and negative. In the Breit frame (and for 
positively charged baryons) , a negative value of the elec- 
tric quadrupole moment can be interpreted as an oblate 



distribution of electric charge. This result agrees with 
lattice estimations, albeit in this case lattice gives very 
noisy results and only for relatively high Q- values. 

As expected, we obtain similar results for the f2~, al- 
though with a different sign coming from the £1 charge. 
The electric quadrupole form factor is non- vanishing and 
negative, and therefore the charge distribution in this 
case also features an oblate shape. This result agrees as 
well with the available lattice data. 



D. Magnetic octupole form factor 

Similar to the electric quadrupole moment, in the Breit 
frame the magnetic octupole moment measures the devi- 
ation from sphericity of the magnetic dipolc distribution. 

In the case of the magnetic octupole, we have to face 
an 1/Q 6 factor when extracting the form factor from the 
electromagnetic current. This entails that, with our cur- 
rent accuracy, we cannot give a reliable value for Gm3(0), 
as is clearly seen in the bottom-right panels of Figs. 5 
and 6. However, in both cases and for both the MT 
and the SD models, we can unambiguously say that the 
magnetic octupole moment is non-vanishing but small, 
and positive (negative for the f2~). We therefore predict 
a prolate distribution of the magnetic dipolc. Unfortu- 
nately, for the magnetic octupole form factor there are 
no reliable lattice calculations to compare with, although 
a quenched calculation [17] suggests a negative sign for 
the A + , in contradiction to our findings. It is very well 
possible that a more elaborate truncation would change 
the sign of our results. However, it is for us very difficult 
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to estimate a priori how the inclusion of, for instance, a 
pion-, resp., kaon-cloud would modify them. 



V. SUMMARY 

We have shown the calculation of the electromagnetic 
form factors of the A and il baryons in the Poincare- 
covariant BSE and DSE framework. This framework has 
as a goal to provide an unified and systematically improv- 
able approach to hadron physics from continuum QCD. 
The calculation presented here uses the Rainbow-Ladder 
truncation of the complete interaction kernel and within 
this truncation scheme we solved self-consistently for all 
the elements in the equations, namely the full quark prop- 
agator and quark-photon vertex. We have performed the 
calculations using two different models for the dressings 
required in the RL truncation, as an attempt to provide 
results which are qualitatively model- independent. 

Our results at u/d quark mass show good agreement 
with lattice calculations and are compatible with the few 
experimental data available for the A. We obtain a neg- 
ative value of the electric quadrupole moment, indicating 
an oblate charge distribution. The sign of the magnetic 
octupole moment is, however, positive, which would cor- 
respond to a prolate magnetic dipolc distribution. In the 
absence of a proper treatment of the current-quark mass 
dependence of mesonic effects or the quark-gluon interac- 
tion in our calculations, we find a weak dependence of the 
electromagnetic properties on the current-quark mass. It 
is, therefore, reasonable that we observe discrepancies be- 
tween our results and lattice calculations for the VL form 
factors. 

This calculation, and especially the magnetic octupole 
form factor, is very sensitive to numerical artefacts and 
for this reason the inaccuracy of the results is sometimes 
significant. Improvements on our algorithms and the em- 
ployment of more elaborate interaction kernels are thus 
desirable in order to verify, in particular, the sign of the 
magnetic octupole moment. 
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Appendix A: Extraction of the form factors 

In Section II we derived an expression for the electro- 
magnetic current in terms of the photon interaction with 
the quarks forming a baryon. On the other hand, the 
form of the current is constrained by Lorentz invariance 
and current conservation to be a linear combination of 
a finite numbers of Lorentz covariants with scalar coeffi- 
cients. These coefficients are the form factors. 

The electromagnetic current for a spin- 3 /2 particle is 
characterized by four form factors P;(Q 2 ) [57, 58]. Its 
expression reads 



J^ Q/3 (P,Q) = P aa ' (P f ) 



(F 1 +F 2 )i 1 ^-F 2 — )S a ? 



(F 3 + P 4 )i 7 " - P 4 



M J 4M 2 
P^(J>-) (Al) 



where P is the Rarita-Schwinger projector 



A+(P) =-(! + . 



r(p) = A + (p)(3r-3^j 



(A2) 
(A3) 



with 7y = Tp 7", Tp the transverse projector (8) and 
the hat denotes a unit vector. Pi and Pf are the initial 
and final baryon total momenta, respectively, Q = P/— P; 
is the photon momentum, M is the baryon mass and 
P = (Pf +Pj)/2. The form factors that are measured ex- 
perimentally are the electric monopole (Ge (Q 2 )), mag- 
netic dipole (Gm 1 {Q 2 )), electric quadrupole (Ge 2 (Q 2 )) 
and magnetic octupole (Gm 3 (Q 2 )) form factors. They 
are related to the P/s via [57] 



2t 

t 



G Bo = 1 + — I (Pi - rP 2 ) - -(1 + r) (F 3 - rP 4 ) ; 



(A4) 



G Ml = (l + y) (Pi + P2) - y (1 + r) (P 3 + P 4 ) , 

(A5) 

G E2 = (Pi - rP 2 ) - \ (1 + r) (P 3 - rP 4 ) , (A6) 
G Ah = (Pi + F 2 ) - X - (1 + r) (P 3 + P 4 ) , (A7) 



with r = Q 2 /AM 2 . It is shown in [57] that if charge and 
magnetic dipole distribution in the baryon is spherically 
symmetric then Ge 2 and Gm 3 must vanish, respectively; 
therefore they measure the deformation of the object. At 
Q 2 = the form factors define the electric charge (e 3 / 2 ), 
magnetic dipole moment (/i3/ 2 ), electric quadrupole mo- 
ment (Q3/2) and magnetic octupole moment (0.y 2 ) of a 
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spin- 3 /2 particle, 



where 



e3/ 2 


= G Ea W , 


(A8) 


M3/2 


= ^(0), 


(A9) 


Q3/2 




(A10) 


o 3/2 




(All) 



Once the electromagnetic current is calculated from 
(24), the form factors can be extracted using the expres- 
sions [58] 



G Eq 

Ge 2 
Gm 3 



s 2 - 2si 



9i 

4o7 (S4 - 2S3) 



8ir 2 VTTT 
16 T 3 



2S3 T 









2si (t- 












-9- 


rs 4 


• 







(A12) 
(A13) 
(A14) 

(A15) 



si = Tr{j^P^P Q P' 3 } 
s 2 = Tr I jm,"«pm | j 

,s 3 = Tr{j^ Q/3 7 £P Q P /3 } 
s 4 = Tr{J"' aa ^} . 



(A16) 

(A17) 

(A18) 
(A19) 
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